Classical dynamics near the triple collision in a three-body Coulomb problem 
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We investigate the classical motion of three charged particles with both attractive and repulsive 
interaction. The triple collision is a main source of chaos in such three body Coulomb problems. 
By employing the McGehee scaling technique, we analyse here for the first time in detail the three- 
body dynamics near the triple collision in 3 degrees of freedom. We reveal surprisingly simple 
dynamical patterns in large parts of the chaotic phase space. The underlying degree of order in the 
form of approximate Markov partitions may help in understanding the global structures observed 
in quantum spectra of two-electron atoms. 
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Our understanding of the overall dynamics of gravi- 
tational or Coulomb three-body problems (TBP) is still 
very limited. This can be attributed to the large di- 
mensionality of the systems, the long range interac- 
tions and the complexity of the dynamics near the non- 
regularisable triple collision pj. Contrary to gravita- 
tional TBP's, where the triple collision is of minor prac- 
tical importance due to the overall stability of these sys- 
tems, the dynamics in three-body Coulomb problems as 
for example two-electron atoms is dominated by the triple 
collision. Wannier's threshold law as well as experi- 
mental results on double ionisation due to single [3j or 
multiple Q| photon processes show electron-electron cor- 
relation effects which are linked to near triple collision 
events 0, |(J . The excitation of two single-electron wave 
packets in two electron atoms has recently been achieved 
experimentally making it possible to study classical col- 
lision events quantum mechanically 0- The importance 
of the triple collision on quantum spectra of two elec- 
tron atoms is also evident from semiclassical treatments 
restricted to collinear subspaces of the full dynamics §. 

We present here for the first time a comprehensive 
analysis of the classical dynamics near the triple colli- 
sion in two electron atoms in the full L = phase space. 
Previous work has been restricted to low dimensional in- 
variant subspaces as for example collinear configurations 

or the so-called Wannier ridge HHH[lJ[l2 We start 
by describing the dynamics at the triple collision in 3 
degrees of freedom. The topology of the space and the 
particle exchange symmetry strongly affect the scattering 
signal for energies below the three-particle break-up en- 
ergy. We derive new scaling laws similar to the Wannier 
threshold law and show for the first time that large 
parts of the phase space is highly structured in terms of 
approximate Markov - partitions. This is an important 
step towards analysing electron-electron correlation ef- 
fects found in quantum properties of two-electron atoms 
Q in terms of the overall phase space dynamics. 

Working in the infinite nucleus mass approximation 
and restricting ourselves to angular momentum L = 0, 
the dynamics may be described in terms of the electron- 
nucleus distances r\ and and the inter-electronic angle 



9. Introducing the scaling parameter, R = (r\ + 7" 2 ) 1 / 2 
and the hyper-angle a = tan _1 (r2/Vi), we employ the 
McGehee transformation 0, 0] to obtain the equations 
of motion for two-electron systems 
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V = — Z(l/ cos a + 1/ sin a) + 1/vT — sin 2a cos 9, 

and Z is the nucleus charge. The scaled differential 
equations are independent of R, and R(t) and H (t) 
are given implicitly through J5J. The triple collision it- 
self has thus been removed from the dynamics and the 
only singularities remaining are the two-body collisions 
at a — 0, 7r/2 which can be regularised by for example 
considering the transformation p a — p a sin 2a. 
There are two fixed points of the dynamics Q), 
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0, p R = ±P , 



where P = [V2(4Z - 1)] 1/2 . These fixed points corre- 
spond to trajectories in the unsealed phase space where 
both electrons fall into the nucleus symmetrically along 
the collinear axis, that is, the triple collision point (TCP) 
with pr = —Pq and its time reversed partner, the trajec- 
tory of symmetric double escape, the double escape point 
(DEP) with pr = Pq. A true reduction in the dimen- 
sionality (from 5 to 4) of the problem is achieved for the 
special initial condition H — for which H becomes a 
constant of motion. The H — - subspace which con- 
tains the two fixed points corresponds to R = or E = 0, 
that is, to the dynamics at the triple collision which is 
equivalent to the dynamics at E = 0! We will come back 
to this remarkable fact when studying the transition from 
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E = to the dynamics for E < 0. There are three in- 
variant subspaces of the dynamics: the collinear spaces 
9 = 7r, pe — (the eZe configuration) and 9 — 0, pe = 
(the Zee configuration) as well as the so-called Wannier 
ridge (WR) of symmetric electron dynamics, a = 7r/4, 
p a = 0. The e-Ze configuration and the Wannier ridge 
are connected at the fixed points. 

Linearising the dynamics near the DEP reveals that 
two of the five eigendirections, v' 1,2 \ lie in the WR with 
eigenvalues A^ 1,2) = -[1 ± yJ{AZ - 9)/(4Z - 1) ] P /4; 
the other three eigenvectors, v^ 3 ' 4,5 ) , lie in the eZe-space 
with A</ 4) = -[1 ± ^/(lOOZ- 9)/(4Z- l)]P /4 and 
A^ } = P , where the subscript D refers to DEP 2]. The 

unstable manifold uffi along the v^ 5 ) direction is ori- 
ented along — pr leading out of the E = subspace; the 
unstable manifold lies in E = and corresponds 
to a trajectory leaving the DEP towards pr = oo. The 
3D stable manifold including S*^ associated with 
(i = 1,2,3) is embedded in the 4D E = subspace. In 
what follows the structure of the Sd is of special impor- 
tance being responsible for the striking features in the 
classical electron - impact scattering signal for total en- 
ergy E < 0. (From time reversal symmetry, one obtains 



E D C B 



for the TCP, A 



-A^ } for i 



1 — 5. The manifolds 



U£' 2 > 3) and 4 4 ' 5) can be obtained from the correspond- 



ing Sp' 2 ' 3 ^ and Up'°' respectively.) 

We analyse first the topology and the dynamics in the 
triple-collision phase space E — with a special focus 
on the structure of the Sd- The topology of the space 
E = is most conveniently studied by considering the 
3D Poincare surface of section (PSOS) 9 = tt in a - p a 
- pa coordinates. From J2J), one obtains that the eZe - 
configuration forms the boundary of the PSOS, see Fig.Q] 
A typical scattering event follows a trajectory coming 
from infinity with pr = — oo, a = or7r/2 and one of 
the two electrons leaving towards infinity with pr — > oo, 
a — ► 7r/2 or 0. After switching to regularised coordinates 
p a , the topology of the eZe - phase space takes on the 
form of a sphere with four points taken to \pn\ = oo, see 
Fig. n . The two fixed points are located at the saddles 
between the arms stretching in forward and backward 
direction along the pr - axis. The WR space connects 
the TCP and DEP along the axis a = tt/4, p a = 0, see 
Fig. 0] it is a compact manifold with the topology of a 
S 2 - sphere where the fixed points form opposite poles. 

For H = 0, we have pn > 0, which leads to a relatively 
simple overall dynamics. Its important features can be 
characterised by the behaviour of the stable/unstable 
manifolds of the fixed points. Various cuts of the Sd in 
the PSOS for fixed p R < P are shown in Fig.ffl The S D 
is for — P < pr < Po bounded by the ID stable manifold 
S^ in the eZe - space, and the 2D WR. That is, the 
Sd is the stable manifold of the WR. What makes the 
evolution of this manifold remarkable, is its behaviour 
at the TCP at pn = — Po, where the phase space itself 
splits into two distinct parts. Starting at pr = Pq and 
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FIG. 1: The PSOS = n of the E = manifold in a-p a -p R 
coordinates. The eZe space forms the boundary of the PSOS, 
the WR connects the TCP and DEP fixed points; below, var- 
ious cuts of the PSOS at fixed pa - values are shown together 
with the Sd- The two arms of the Sd stretching from the 
WR towards the S^ on the eZe boundary are shown as full 
and dashed line, respectively. (The Sd in C and D is drawn 
schematically to enhance important features.) 



going towards decreasing pr - values which corresponds 
essentially to an evolution of the Sd backwards in 
time, the Sd undergoes the usual stretching and folding 
mechanism typical for an unstable manifold in compact 
domains. The stretching and folding is here facilitated 
by an overall rotation of the space around the WR - 
axis a — 7r/4, p a — and a certain "stickiness" near 
the binary collision points a = or 7r/2, see the cuts B 
and C in Fig. ^ UM- PR ~^ ~Po, the phase space 
develops a bottleneck whereas the Sd stretches over the 
whole phase space 5 times by now. As pu decreases 
further passing through — Po, the Sd is cut at the TCP 
into distinct parts, see D in Fig. 2] Points close to the 
TCP will leave the fixed point along the stable manifold 



S^ in eZe space towards pr 



In each arm 
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exactly 5 pieces of the Sd are connected to the Sj, for 
Pr < — Po forming two loops and one connection to the 
eZe boundary. The S^ itself is thus also a boundary 
of the S d without being part of it. There are therefore 
two main routes towards the DEP for electrons coming 
in from pr = — oo close to the eZe - boundary: the first 
route approaches the TCP on one of the 5 leaves of the 
Sd close to the stable manifold S^ and moves then 
along the WR towards the DEP; secondly, trajectories 
can approach the DEP 'directly' by moving on the Sd in 
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FIG. 2: The time-delay signal for 9 = tt/2 and £i = 0.2; 
five distinct peaks appear in the 'dips' associated with close 
encounters with the triple collision point. Note that delay 
times are shown on logarithmic scales. 



the vicinity of the . This twofold approach is essen- 
tial for understanding the dynamics in the E < - space. 

We now turn to the full dynamics in the 5D E < 
phase space. In the following we analyse one-electron 
scattering signals where electron 1, say, starts at r\ = oo 
with energy E\ and fixed total energy E = — 1. We record 
the time delay of the outgoing electron. (Note that fixing 
E = —1 is sufficient as changing the total energy amounts 
to a simple scaling transformation in the dynamics [To|). 
A smooth transition of the dynamics from E < towards 
E = is achieved by considering the limit E\ oo, 
E2 — ► — 00 and E1+E2 = — 1. The inner electron is then 
initially bounded infinitely deep in the Coulomb well and 
interactions between the incoming and bound electron 
will take place at R — > and thus H — > 0. The dynamics 
in E — is in this sense equivalent to the dynamics at 
the triple collision point R = and the phase space flow 
at E = can be smoothly continued to the flow in the 
full E 7^ space. Trajectories close to H = will follow 
the dynamics in the triple collision space except near the 
fixed points where the flow close to the manifold E = 
is perpendicular to that manifold along the v' 5 ) = —pr 
direction. 

In Fig. |21 a typical signal for the scattering time r 
is shown, here for 8 = ir/2 and E\ = 0.2, where 9 is 
the inter-electronic angle for n — > 00. (We have cho- 
sen zero angular momentum of both electrons initially). 
ip 6 [— 7r,7r] parametrises the starting point along one 
revolution of the trajectory of the inner electron. Disre- 
garding the finer structures of the signal at the moment, 
two main features emerge, a dip in the scattering time 
around <p = —1.0 and a chaotic scattering interval (CSI) 
around 0.2 < cp < 2.2. The dip can be associated with 



trajectories close to St near which move towards the 
TCP. Most of these events will lead to immediate, fast 
ionisation of one of the two electrons due to the large 
momentum transfer possible near the triple collision sin- 
gularity causing the dip in the scattering time. Note 
that the 2D stable manifold St is completely embedded 
in the eZe configuration; it is not possible to reach the 
TCP from outside the eZe region. 

Chaotic scattering occurs, on the other hand, if trajec- 
tories approach the DEP close to the stable manifold Sd 
via the direct route near S^ 3) . The 3D stable manifold 
of the DEP is part of the E — space and the DEP can 

thus be reached for E < only in the limit E\ — > 00. 

(3) 

Orbits close to the S D will, however, also come close to 

the DEP where they either follow Up leading to ioni- 
cs) 

sation or follow U D along the pr axis perpendicular to 
the H = - manifold. In the latter case, pn changes sign 
and electron trajectories fall back towards the nucleus 
having equidistributed their momenta in such a way that 
chaotic scattering becomes possible. The DEP thus acts 
as a turnstile for entering into a chaotic scattering region 
of the full phase space. By time-reversal symmetry, the 
TCP acts as the exit gate for single electron ionisation. 

A closer analysis of the strongly fluctuating signal in 
the CSI reveals a series of dips flanked by singularities in 
the delay time on either sides, see Fig. [21 The dips cor- 
respond to orbits which after having spent some time in 
chaotic motion leave the chaotic region by coming close 
to the TCP along the stable manifold Srp . The borders 
of these intervals are given by orbits escaping asymptot- 
ically with zero kinetic energy of the outgoing electron. 
All these escaping regions can be labelled uniquely by 
a finite binary code reflecting the order in which binary 
near-collision events take place after entering and before 
escaping the chaotic region. The mechanisms leading to 
escape and to the existence of a binary symbolic dynam- 
ics forming a complete Smale horse-shoe are well docu- 
mented and understood for the eZe configuration |9Ul0|. 
It is surprising that this mechanism holds in principle 
also far away from the eZe regime down to initial inter- 
electronic angles of the order 9 m 7r/4. 

There is one major difference between the scattering 
signal for the collinear eZe space and the dynamics for 
9 7^ 7r; additional structures emerge at the centre of the 
dips, see Fig. [2J When enlarging the regions near the 
dips both at the primary dip and in the CSI, one finds 
five separate peaks and, on further magnification, each 
of these peaks breaks up into a chaotic scattering signal 
similar to the primary pattern. One finds in this way a 
whole sequence of self-similar structures where dips give 
birth to CSI's which in turn have 2nd generation dips 
containing 5 peaks etc. Such a sequence leading up to 
the second generation CSI is shown in Fig. [21 

These phenomena can be explained in terms of the dy- 
namics in the triple collision space H — discussed ear- 
lier. The dips represent close encounters with the TCP; 
whereas in the eZe - space, the only route away from 
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FIG. 3: Energy scaling of both the width of the primary CSI 
Acs/ (see Fig.[5J for different 6 values as well as the width of 
the 5 peaks As p (here for 6 = 7r/2) is shown; the scaling agrees 
with n = 1.0558932... for Z = 2. To the right, the width, 
Acp, of the centre peak of the 5 peaks is shown as a function 
of n — 8 (for fixed E\) with scaling exponent v — 6.223573.... 



(3) 

the TCP is along the unstable manifold Uj, , for ^ 7r 
another route opens up: escape along the WR and thus 
along parts of the Sd leading back to the DEP and to re- 
injection into the chaotic scattering region! The existence 
of the 5 peaks in each of the dips is thus a clear manifesta- 
tion of the stretching, folding and cutting mechanism of 
the Sd in the H — dynamics. A closer analysis reveals 
that the centre peak (CP) of the 5 peaks corresponds to 
the part of the S d connected to the WR; the outer peaks 
contain orbits which move away from the WR after pass- 
ing the TCP and before entering the DEP region. By 
repeatedly moving from the DEP to the exit channel, the 
TCP, and then along the WR back to the DEP, it is pos- 
sible to create increasingly longer cycles of chaotic scat- 
tering events. We do thus encounter here a rather curious 
dynamical feature, namely a Smale-horseshoe, whose en- 



trance and exit points are short-circuited by a degenerate 
heteroclinic manifold, the WR, connecting the two fixed 
points sitting at these turnstile-gates. 

The features described above imply scaling laws for 
the width A of the CSI's for E\ — > oo. In this limit, 
trajectories which will enter the chaotic scattering region 

(5) 

along U D need to come closer and closer to the DEP. The 
size of this phase space region near the DEP is limited 
by ejection along t/p . A simple estimate yields A ~ 

E'i ^ with /i = A^/A^J . The universality of this scaling 
law for all CSI's is shown in Fig. [3] and demonstrates 
that the DEP is indeed the sole entrance gate into the 
chaotic scattering region. The exponent \x is the so-called 
Wannier exponent which plays a crucial role in Wannier's 
threshold law for E > Q and is also of importance for 
quantum resonance widths [Tfij . A similar scaling law can 
be deduced for the CP; the phase space volume which can 
be transfered from the TCP to the DEP along the WR 
scales in the eZe - limit 9 — > it like Aqp ~ (vr — 6) v E^ 



with 



A 



(3) 



/Re[X^] = X^'/RelX^}, see Fig. ED 



v(3) 



Our analysis suggests that the Markov partition lead- 
ing to the symbolic dynamics of the eZe - space remains 
largely intact over a substantial part of the L — phase 
space with the only modification that the number of sym- 
bols increases from 2 to 2 + 5 = 7. The existence of such a 
partition in this high-dimensional problem is not obvious 
and may be the key for explaining approximate quantum 
numbers observed in two electron atoms in terms of the 
classical dynamics @- 
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